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! 1 Statement and discussion of obtained results 

. Any value x e M is representable as a continued fraction (CF) 
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, where ap G Z (if oq — 0, then we omit it), e N for alH > 1. Let i be a fixed position in the CF expansion 

^ . and let Pi{k) stand for the probability that = k with x randomly chosen in the segment [0, 1). Hereinafter 

we understand the random choice as a realization of the uniform distribution. It is well known that according 
to the Gauss-Kuzmin theorem, 

hmP.(fc)=ln(l + ^^)/ln2. (1) 

The roots of the quadratic equation x"^ + px = q, where p,q E I1, A = + 4g > 0, are expandable in a 
periodic CF. Let x+{p^ q) = \/A/2 — p/2. The fractional part of this value {a;+(p, q)} = x+{p^ q) — [x+(p, q)\ 
has no pre-period in the CF cxpansiorQ: {x+{p,q)} = [[oi, . . . ,aT(p.g)]]. Here T(ji,q) is the length of the 
period of the CF for x+{p, q); denote the period itself (ai, . . . , aT{p,q)) by a{p, q). 

Let us adduce the Arnold conjecture about the Gauss-Kuzmin statistics for periods of CF of quadratic 
irrationalities ([H problem 1993-llB], mHH): 
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Conjecture 1. Choose an integer point {p, q) in the circle of radius R with the probability proportional 
to T{p,q). Further, let us randomly choose Oi from the corresponding set a{p,q) (V.I. Arnold performs this 
two-step procedure immediately by the random choice from the union of sets a(ji, q) ). Then the probability 
that Oi — k tends to ln(l + ■^jrpsy)/ as R ^ oo. 

Unfortunately, for some reasons, the proof of this conjecture performed by V.A. Bykovskii and his fol- 
lowers was not published. Below we describe the simple proof of the following theorem: 

Theorem 1. Let w G (0, 1) be a fixed value. Let us randomly choose a fixed integer point {p, q) in the 
circle of radius R. Further, let us choose from the corresponding set a{p, q) the value Oi with the probability 
proportional to w*. Then the probability that at — k tends to (1 — w) X^i^i Pi{k)w'''~^ as i? — > oo. The limit 
of this sum with w ^ 1 is the Gauss-Kuzmin statistics ln(l + -^^^)/hi2. 

The latter part of this assertion, evidently, follows from the Kuzmin theorem ([1]) and the regularity of 
the Abel summation method [5] . We did not succeed to prove a variant of Theorem 1 with w initially equal 
to 1. 

Note that with the random choice of (p, q) in the circle of radius R the number x+ (p, q) is not necessarily 
an irrational real value. In accordance with Theorem 1 the probability of this event tends to zero as i? — > oo. 
Let VLr stand for the set of pairs (p, q) for which this event does not take place with fixed R. 

In papers [3l |4] V.I. Arnold experimentally studied the value 

f(i?)= T{p,q)/\nj,\. 

He established that 

f{R) - const R asR-^oo. (2) 

Let To{q) = T(0, q) be the period of the CF for the square root of q. Put fo(g) = Y.q=i To{q)/Q. Some 
experts in the number theory studied the mean value of the mentioned kind experimentally; this enabled 
them to state the following conjecture: Tq{Q) ~ const-^Q ln^"((5), where a > (see ^ and references 
therein) . 

Note that one can easily upper estimate To(Q) [Z]: 

fo(Q) < const^. 

In [S] (see also [B]) E.P. Golubeva proves that the left-hand side of this inequality is asymptotically small in 
comparison with the right-hand one; moreover, if the extended Riemann conjecture is true, then the order 
of their difference is not less than ln((5)'"^^^. The bound 

fo(g) > const v/Qln""(Q) (3) 

would allow us to answer the famous Gauss question [9j about the growth of the number of classes of real 
quadratic fields. The known lower bound for Tq{Q) is far from the right-hand side of (0]). In [6] E.P. Golubeva 
proves only that Tq{Q) > const ln(Q). 

Problem 5 listed in [lOj implies the estimation of the growth rate of elements of the period of a CF. Let 
us adduce its exact statement [lOl page 7]. 
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Let a{p,q) = Eaiea(p,g) "'/^(P' ^(-R) = E(p,g)gOfl "The problem is to evaluate the 

growth rate of A{R): is A greater than Ci?" for some positive C, a? Or is it smaller than some C(ln R)" ?" . 
Let us prove the following assertion. 

Theorem 2. Let us choose ai in the same way as in the statement of Conjecture 1; let A'p, stand for the 
mean value of the chosen number, that is, 



E(p.g)60„ T{p, q) 



// correlation 0) or at least that T{R) > const R (In R) ^ analogous to inequality (0) takes place, then we 
have 

A'jj < C(lni?)". 



The main idea of the proof of Theorem 1 is the representation of the probability under consideration 
as the Riemann integral. By applying the Weyl theorenj^ to the uniformly distributed sequence a;+ {p, q) in 
Section 2 we easily prove the desired assertion. 

The idea of the application of the Riemann integral sums was used implicitly for the proof of a similar 
correlation in the case of rational values with fixed denominator [12, Theorem 4.5.3.E]. However, in this 
case one usually applies another technique. An analog of Conjecture 1 for rational values p/q was proved 
by M.O. Avdeeva and V.A. Bykovskii in [T3]. The asymptotic validity of the Gauss-Kuzmin statistics in 
the case of a fixed denominator follows from results obtained by H. Heilbronn and J.W. Porter (see [14]). In 
[T5] A.V. Ustinov establishes the limit statistics for finite fractions whose divisors and denominators belong 
to an arbitrary expanding domain. 

The proof of Theorem 2 is based on the upper estimate for X)a ea(p q) ^« ^® ^ certain simple function of the 
discriminant A. We prove this inequality (possibly, known by experts) with the help of a simple geometric 
construction that represents a symbiosis of the "Nose-Hoover" algorithm [2] and an explicit technique for 
defining values of an integer- valued sign-indefinite quadratic form [16j . 

Note that sign-indefinite quadratic forms correspond to hyperbolic operators. Symmetric properties of 
hyperbolic operators were studied in geometric terms by F. Aicardi ([17], see also [T51[Tn]). The algebraic 
approach to this topic is described in well-known papers [201 [H] . 

This paper is written under the bright impression of lectures delivered by V.I. Arnold [521 1231 121] and a 
result of the performed numerical experiments. 



2 Proof of Theorem 1. 

Let b be an arbitrary number from the segment [0, 1]. Recall (05]) that a sequence is said to be uniformly 
distributed, if among its first N elements the quotient of terms whose fractional part is less than b tends to 
& as ^ oo. 

^The reference to the Weyl theorem was done by a reviewer of the journal "Functional Analysis and Its Application", the 
proof of the first variant of Theorem 1 adduced in is more awkward. 
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Figure 1: Domains considered in the proof of Lemma 2. The rays that bound the domain 5" are tangent hues 
to the paraboHc curve q = — Indicated are sets of points M2o,i,i/2 H ^ ^20,1,1/2 -^20 1 1/2- 

Lemma 1. Let us enumerate numbers x+{p, q), (p, q) E 1? , + Aq > Q in ascending order of the distance 
from the point {p, q) to the origin of coordinates (in the case of equal distances the numeration is arbitrary). 
The obtained sequence is uniformly distributed. 

Since x^{p,q) — p, g) G Z, in order to prove the lemma, suffice it to consider the terms of the 
sequence that correspond to points in the right half-plane. Lemma 1 easily follows from Lemma 2. 

Lemma 2. Let n is a fixed integer number. The subsequence of the sequence from Lemma 1 that corresponds 
to points {p, q) such that n<q/p<n + l,p>0 is uniformly distributed. 

Proof of Lemma 2 Conditions n < q/p < n + 6, p > 0, 6 G [0, 1] define the sector Sn,n+b hi the plane 
{p,q). Instead of it, let us consider another domain S!^ ^^i^. The domain 5^ is defined by inequalities 
np + n"^ < q < {n + b)p + (n + fe)^, with p > — 2n and inequalities — < q < {n + b)p + (n + 6)^ with 
—2n — b < p < -~2n (see Fig. [T|). One can easily make sure that the rays that define the boundaries of the 
domain S' are equiscalar lines of the function a;+(p, (7). Evidently, the condition {x+(p, (7)} < b for points 
from S'„ „^i is equivalent to their belonging to the subsector 5*'^^ 

Let MR,n,b {M'j^ „ J,) stand for the set of integer points located inside the sector Sn,n+b (the domain 
^'n7i+b) ^^'^ the circle of radius R centered at the origin of coordinates. Let Mii^n,b = Ma^n,! \ MR^n,b, 
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Evidently, for i? ^ cx) we have |M;^„_j, n Mfl,„,b|/|Mfl„„,i| ^ b, \M'ij^^^^f, n Affl,„,b|/|Mfl^„,i| ^ 1 - fe, 
whence we obtain the assertion of Lemma 2. □ 

The assertion of Lemma 1 follows from the inequality | Un |n|<7v -^^i?,n,i|/|-^'^fl,+ | > 1 — 1/-^, where M^^^ 
is the collection of all integer points located in the half-circle of radius R. 

Let us now immediately prove Theorem 1. Let i,k G N, let /^(s) stand for the indicator function 
of the set of points x from the segment [0, 1] for which the i-th position in the CF is occupied by the 
number k. This set and its complement are representable as unions of countable numbers of intervals [26j . 
Consequently, the function I^x) is countably continuous and therefore it is Riemann integrable. The function 
fk,w{x) = (1 — lu) X^i^i -^k{^)'^^~^ the same property. The integral of this function 

„1 oo 

/ (l-w) VP,(fcK-i (5) 

equals the probability of the event A: aj — k with the random choice of the number x in the segment [0, 1). 
Here j is chosen randomly, namely, Prob(j = i) = (1 — i — 1, 2, 3, . . .. 

We can calculate the Riemann integral ([5]) with the help of the Weyl theorem [5S]. In accordance with 
this theorem the frequency of the event A for the first N terms of a uniformly distributed sequence tends 
to its probability. Taking into account the fact that the CF for {a;-|_(p, q)} has no pre-period, we obtain the 
assertion of Theorem 1. □ 

3 Proof of Theorem 2 

We remind [2] that the "Nose-Hoover" algorithm for finding a CF of a real number x, x > 0, is reduced to 
the geometric method which constructs the boundary of the convex shell of the set of integer nonnegative 
points {u,v) located above (below) the straight line v — xu. Let Bq stand for the vector (0, 1), ei = (1,0). 
Put 

e„+i = e„_i + a„_ie„, 71 = 1,2,..., (6) 

where a„_i is the maximal integer such that the vector e„+i lies below (for odd n) or above (for even n) 
the straight line v = xu. Klein noted that a„, n = 0, 1, . . ., coincide with partial quotients of the CF of the 
number x. The geometric algorithm results in two polylines which represent parts of sails of the CF. 

Let x > \. We need a slightly modified algorithm which results in one infinite polyline L, originating 
at the point (1, 1), whose segments are the vectors a„_ie„, rt = 1, 2, . . .. Evidently, e„+i + Cn = e„ + Cn+i- 
Consequently, if in the standard "Nose-Hoover" algorithm ^ we add the extra vector e„ and thus go out 
of the line v = xu (i. e., we consider the vector e„_i -I- (a„_i -t- l)e„), then we get the first integer point on 
the segment of the polyline of another sail constructed at step n + 1. This fact justifies a simple geometric 
algorithm which constructs the polyline L. 

We begin the construction process with the point (1, 1); it is convenient to connect it with the origin of 
coordinates by the segment which does not enter in L. As the "constructive" term at the first (nth) step 
we choose the vector Ci (the vector e„). We add this vector till we go out of the line v = xu. The newly 
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Figure 2: The gradual "Nose-Hoover" algorithm for x = \/2 



added vector e„+i is directed from the origin of coordinates to the point obtained as a result of the latter 
addition up to the step out of the line. See Fig.[2]for the first segments of the polyline for x = \/2. The zero 
partial quotient (i. e., [x\) differs from zero. In order to provide the correspondence to the indices of partial 
quotients, it is also convenient to begin the numeration of segments of the polyline L with zero. 

Let us consider the approximation of a number x by mediants. Let /i = ^7, /2 = f| be irreducible 
fractions such that wi,ui,U2,U2 > and /i < /2. Let us define the operation I of finding the mediant 
(the "insertion" operation) by the formula /i 1/2 = ^, where v — vi + V2, u ^ ui + U2- Earlier we used 
the denotation introduced by A. A. Kirillov. Evidently, in the geometric representation of the fraction ^ as 
the integer vector (m, w) the operation J, corresponds to the addition of vectors {ui,vi) and {u2,V2)- The 
obtained diagonal of the parallelogram appears to be "inserted" between its sides. 

Let X G (/i, /2), /a = fi i /2- Since /s G (/i, /2), one can treat /a as an approximation of the number x. 
If a; € (/i, /s), then we put /4 = /i i /a, otherwise we do /4 = /s i /2- The process of the approximation of 
the number x by mediants consists in the repetition of these operations for the corresponding intervals. 

The condition a; > means that x £ (f; 5)- The initial approximation of x for such an interval is the 
fraction j. One can easily see that for the mapping 

y 

the fraction — the point (u, v) 
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the geometric representation of the algorithm for approximating the number x by mediants is completely 
identical to the algorithm for constructing the polyline L. The technical distinction consists in the following 
fact: earlier we constructed each segment of the polyline L "at once", but now it "grows gradually" due to 
the stepwise addition of the next vector e^. 

Further we consider a simplified version of this algorithm for the quadratic irrationality x^(j)^ q). We will 
need it for the proof of the next assertion. 

Lemma 3. For n G N put D(n) — X^i^^ '''i'^ ~ u"^): where [-J is the integer part, T{m) is the quantity of 
divisors of the number m. We have the inequality 

^ a, < /(A/4), where A = p^ + Aq, (7) 

ai£a{p.q) 

f{n) ~ 2D{n) + T(n), if n is integer; (8) 
f{n) — 2 r(ri — otherwise. (9) 

iG{l,3,...} 

<4n 

In the case of an odd period T{p, q) one can improve this bound, namely, divide the right-hand side of 
inequality ^ by two. 

Proof of Lemma 3 Note that {x^[p, g)} = {x^[p + 2, q — p — 1)}, therefore, without loss of generality, we 
assume that p — or p = 1, here x^{p, q) > I (the special case a;+(l, 1) is evident) and the second root of 
the quadratic equation is negative. 

Let us apply to x+{p, q) the gradual "Nose-Hoover" algorithm described above. Let (u„, «„) stand for the 
nth integer point on the polyline L, n = 0,1,2,.. .. We obtain it at the nth step of the gradual "Nose-Hoover" 
algorithm. In Fig. [2 for example, iuo,vo) = (1,1), ("1,^1) = (1,2), {u2,V2) = (2,3), (u3,W3) = (3,4), 
(m4, ^4) = (5, 7), etc. 

Let us draw one more vector which originates from the point (u„, w„); denote it by e^. The vector is 
defined by the condition + e; + e" — 0, where is the vector which originates from the point (M„,f„) 
and goes along the polyline L, e" = —{un,Vn) (we considered this vector earlier). The collection of vectors 
{e^,ei,e"}, where each one is defined accurate to the multiplier (—1), is a superbasis [16j . This means that 
any pair of these vectors generates the whole integer lattice. The transition from the point (u„, t;„) to that 
{un+i,Vn+i) corresponds to the replacement of one superbasis with another one; the latter differs from the 
initial superbasis only in one of three its elements. Thus, for the polyline L in Fig. [2] the initial superbasis 
is {(0,1), (1,0), (1,1)}, the 1st one is {(1,1), (1,0), (1,2)}, the 2nd one is {(1, 2), (1, 1), (2, 3)}, the 3rd one is 
{(1,1), (2, 3), (3, 4)}, etc. 

Let us extend vectors up to the rays which originate at the points (««,««). As a result, the first 
quadrant appears to be divided onto disjoint connected domains (see Fig. [2]). Let us associate each vector 
iun,Vn) with a domain, whose boundary contains all points of polylines which include the vector (m„,d„) 
in the corresponding superbasis. In the picture these domains look like narrow "crevices" located at the 
north-east of the points (zi„,w„). We associate the domains which border on the coordinate axes with the 
unit vectors of these axes. 
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Figure 3: The arithmetic progression rule 
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Table 1: All values of {a,b,h) for the polyline represented in Fig. [2] 



We have 

{un, Vn) is located above the line v — x+{p, q)u v'^ + pVnUn > quf^- (10) 

Let us write the values of the quadratic form v'^ + pvu — qu^ at the points in the corresponding 

domain. See Fig. [5] for the result obtained in the case p — 0, q — 2 (the upper right corner of the figure 
contains the value calculated at the point (5, 7)). From (fTO)) it follows that the points of the polyline L belong 
to boundaries of domains with positive and negative values. The collection of superbases which correspond 
to these points is called the river of the quadratic form [16] . 

Let us adduce the main result of the theory of quadratic forms which enables one to make the calculation 
of the period of a CF much easier. Using the values of the quadratic form at three vectors from any 
superbasis, one can easily restore the CF. To put it more precisely, the following correlation [16j is true. Let 
domains associated with values a,b,c,d of the quadratic form be located in accordance with Fig. [31 Then 
values d,a + b,c form an arithmetic progression. Moreover, if h is the common difference of this progression, 
then 

ab^ A/A. (11) 

Therefore, using the values (a, b) for two neighboring domains with different signs separated by the 
polyline L, and the value h, one can unambiguously restore all subsequent and previous values of partial 
quotients of the CF, moving along the river of the quadratic form (see the figures in Chapter 1 of book [16J). 

Let us now immediately prove the inequalities They follow from correlation (|lip. Assume that at 
step ni we get the same parameters of the arithmetic progression as those obtained at step 0. Such a step 
exists because the algorithm is invertible. Let li stand for the number of the polyline segment which contains 
the point (u„j,w„j). Then a part of the CF 01,02, . . . , oj^ becomes periodic. In addition, the number li is 
even, because the points (uo,vo) and (u„j,w„j) are located to one side of the polyline L. In the case of an 
odd period the sequence oi, 02, . . . , aj^ contains at least two repeating subsequences. In accordance with the 
algorithm we have X^'Li o.i = ^i- Therefore, suffice it to estimate the number of all possible triplets (a, b, h) 
which satisfy (fTTj). The obtained bound obeys formulas (|8|9|1 . The multiplier 2 in the right-hand sides of 
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these formulas appears, because we have to take into account various signs of /i, and the term t(A/4) does 
when we consider the case h — 0. Lemma 3 is proved. □ 

Remark In the proof of Lemma 3 we estabhshed that CF for {x+{p, q)} has no pre-period. 

Let us complete the proof of Theorem 2. Formally speaking, the Dirichlet theorem (the equality 
Eu=i ■''('")/"- = li^"- + (27 — 1) + 0{l/y/n)) does not imply that D{n) ~ y^lnn, because the values of 
the function t are not uniform. A more accurate bound is 

D{n) = 0(ln^(n)V^) 

(see [7] and references therein). This, evidently, implies that f{n) — 0(ln'^(n)Y^). Taking the sum over 
all possible values of {p, q), we obtain that the numerator of fraction (jH) equals 0{R'^ ln'^(i?)), whereas the 
denominator by assumption is greater than const (InR)^^. Theorem 2 is proved. □ 
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